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ABSTRACT 

We present a review of the current status of very high energy 7-ray astronomy. The development 
of the atmospheric Cherenkov imaging technique for ground-based 7-ray astronomy has led to a rapid 
growth in the number of observatories. The detection of TeV 7-rays from Active Galactic Nuclei was 
unexpected and is providing new insights into the emission mechanisms in the jets. Next generation 
telescopes are under construction and will increase dramatically the knowledge available at this extreme 
end of the cosmic electromagnetic spectrum. 

Subject headings: gamma rays: observations — galaxies: active — BL Lacertae objects: general — 
supernova remnants — pulsars: general — intergalactic medium 
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1. INTRODUCTION 

Very high energy (VHE) 7-ray astronomy, (defined here 
as observations at energies above 300 GeV and below 
100 TeV) became viable with the development of the at- 
mospheric Cherenkov imaging technique which crossed the 
vital detection threshold with the detection of the Crab 



Nebula ten years ago (Weekes et al. 1989). Although at 
the periphery of the observable electromagnetic spectrum, 
the VHE band must now be considered a legitimate astro- 
nomical discipline with established sources, both steady 
and variable, both galactic and extragalactic, a growing 
number of observatories, and the promise of significant ad- 
vances in detection techniques in the next few years. The 
confirmation of the detection of the Crab Nebula by more 
than eight groups over the past decade has given credi- 
bility to the existence of sources of TeV 7-rays and led 
to a rapid improvement in the sensitivity of ground-based 
7-ray detection techniques. 

The challenge therefore in this short review, which is 
aimed at the general astronomical community, is to con- 
vince the reader that there are viable methods of detect- 
ing 7-rays of energy 300 GeV and above from the ground, 
that a population of credible sources of various classes have 
been detected, and that these detections make a significant 
contribution to the astrophysics of high energy sources. 
Although, as we shall see later, there are plans to extend 
the ground-based detection techniques down to energies of 
20 GeV, and there are viable air shower experiments that 
operate at energies of 50 TeV and above, this review will 
focus on observations in the 300 GeV to 30 TeV energy 
range; this has been the most studied energy band be- 
cauac it is caaily acccaaiblc to the atmoaphcric Cherenkov 
imaging technique. 



2. TECHNIQUES AND INSTRUMENTATION 

2.1. Space telescopes 

The physics involved in the interaction of photons with 
matter is well established. At energies above 10 MeV the 



predominant interaction is pair production in which a 7- 
ray is converted into an electron-positron pair in the pres- 
ence of a nucleus: E 7 — v m e +c 2 + m e -c 2 . The electron 
and positron carry information about the direction, energy 
and polarization of the primary 7-ray and hence detection 
methods focus on the observation of these secondary par- 
ticles. The interaction length of photons in matter is 30 g- 
cm 2 and the earth's atmosphere is 1030 g-cm 2 thick, so 
the preferred (only) method of unambiguously detecting 
high energy 7-rays is from space vehicles (from balloons 
initially but now from satellites). 

The basic elements of a space-borne 7-ray detector are 
therefore (1) a particle detector in which the 7-ray inter- 
acts and the resulting electron pair tracks are recorded; 
(2) a calorimeter in which the electrons are absorbed and 
their total energy recorded; (3) an anti-coincidence shield 
surrounding the tracking detector which registers (and re- 
jects) the incidence of charged cosmic-ray particles (which 
are about 10,000 times more numerous than the 7-rays). 
In most of the telescopes flown to date, the tracking detec- 
tor has been a spark chamber, the calorimeter a Sodium Io- 
dide crystal and the anti-coincidence detector a thin sheet 
of scintillator. The resulting telescope is a sophisticated 
device which can unambiguously identify 7-rays, has an 
energy resolution of about 15%, an angular resolution of 
1° and a field of view of 20° to 40° half angle. Unfor- 
tunately this sophistication does not come cheaply; the 
effective collection area is only a small fraction of the to- 
tal telescope area so that in the largest telescope flown to 
date, the hugely successful Energetic Gamma Ray Exper- 
iment Telescope (EG RET) on the Compton Gamma Ray 
Observat ory (CGRO) ( Thompson et al. 1995| ; Hartman ct 
al. 1999J), the effective collection area was only 1,500 cm 



(about the size of two pages of this journal!) whereas the 
actual physical instrument was about the size of a compact 
car. All of the 7-ray telescopes flown to date (the SAS-II 
telescope in 1973, the COS-B telescope in 1975, EGRET 
in 1991) have had the same functional form; the next gen- 
eration space telescope, the Gamma-ray Large Area Space 



VHE Gamma-Ray Astronomy 



Teles cope (GLAST ), scheduled for launch in 2005 flGehrclsj 
& Mi chclson 1999| ) has the same general features but will 
use solid state detectors with a factor of 10-30 improve- 
ment in sensitivity. 

The Third EGRET catalo g (based on some four years 
of observation by EGRET) (Hartman et al. 1999) con- 
tains a listing of more than 250 sources of >100MeV j- 
rays, more than half of which are unidentified with any 
known astronomical object. In addition to a detailed map 
of the diffuse emission along the galactic plane, there is 
evidence for more than 100 galactic sources, a small num- 
ber of which have been identified with pulsars and possibly 
with supernova remnants. The bulk of the sources away 
from the plane are extragalactic and have been identified 
with active galactic nuclei (AGNs), almost all of which are 
blazars. Some of the sources are variable indicating that 
the high energy 7-ray sky is a dynamic place. Many of 
the identified sources (pulsars, AGNs) have very flat spec- 
tra and have luminosities that peak in the high energy 
region of the spectrum. The EGRET sensitivity extends 
to 10 GeV but is limited by the calorimeter (at the high- 
est energies the cascade from the electron-positron pair is 
not contained and charged particles from the cascade can 
reach the anti-coincidence detector, causing a veto of the 
event ). Since the source flux almost invariably decreases as 



energ y increases, it is only possible to extend observations 
by building larger telescopes. Even for a flat spectrum 
source (power law with differential spectral index -2.0), 
the power sensitivity of EGRET falls off with energy. 

2.2. Ground-based Telescopes 

The earth's atmosphere is as opaque to photons of en- 
ergy > 300 GeV as it is to photons in the CGRO range 
(100 keV to 30 GeV). However, at these higher energies the 
effects of atmospheric absorption are detectable at ground 
level, either as a shower of secondary particles from the re- 
sulting electromagnetic cascade or as a flash of Cherenkov 
light from the passage of these particles through the earth's 
atmosphere. As in the CGRO range, 7-ray observations 
are severely limited by the charged cosmic particle flux 
which gives superficially similar signals at ground-level 
and, for a given photon energy, is 10,000 times more nu- 
merous. It is not possible to veto out the charged cosmic- 
ray background with an anti-coincidence shield. As such, 
it might seem impossible to do 7-ray astronomy with such 
indirect techniques. However there are small, but signifi- 
cant, differences in the cascades resulting from the impact 
of a photon and a proton on the upper atmosphere; also 
the electromagnetic cascade retains the original direction 
of the photon to a high degree, and the spread of secondary 
particles and Cherenkov photons is so large so that a sim- 
ple detector can have an incredible (by space-based 7-ray 
detector standards) collection area. 

In practice the most successful detectors are atmo- 
spheric Cherenkov imaging telescopes (ACITs) which 
record the images of the Cherenkov light flashes and 
which can identify the images of electromagnet ic cascades 
from putative sources with 99.7% efficiency ( Aharonian| 
fc Akcrlof 19971 |Ong 1998|) . Originally proposed in 1977 
( Wcckcs fc Turver 1977| ), the technique was not demon- 
strated until ten years later when the Whipple Observa- 
tory 10-m reflector (Figure fll) was equipped with a primi- 
tive imaging camera and used to detect the Crab Nebula 



( Wcekes et al. 1989| ). The technology was not new (ar- 
rays of fast photomultiplier tubes in the focal plane of 
large optical reflectors with readout through standard fast 
amplifiers, discriminators and analog-to-digital convert- 
ers) but the technique was only fully exploited in the past 
decade. Compared to high energy space telescopes such as 
EGRET, ACITs have large collection areas (>50,000m 2 ) 
and high angular resolution (^0.1°). ACITs also have rea- 
sonably good energy resolution (~ 20 — 40%), but small 
fields of view (FOV) (<5°) and a background of diffuse 
cosmic electrons which produce electromagnetic cascades 
identical to those of 7-rays. ACITs also have low duty 
cycles (<10%) because the Cherenkov signals are faint 
and produced at altitudes of several kilometers, requir- 
ing cloudless, moonless skies for observations. Most of the 
results reported to date have been in the energy range 
300 GeV to 30 TeV 

In recent years VHE 7-ray astronomy has seen two ma- 
jor advances: first, the development of high resolution AC- 
ITs has permitted the efficient rejection of the hadronic 
background, and second, the construction of arrays of AC- 
ITs has improved the measurement of the energy spectra 
from 7-ray sources. The first is exemplified by the Whip- 
ple Observatory 10-m telescope with more modern ver 



sions, CAT, a French telescope in Pyrenees (Barrau et al 



1998), and CANGARO O, a Japanese-Aus tralian telescope 
in Woomera, Australia (Hara et al. 1993). The most sig- 



nificant examples of the second are HEGRA, a five tele- 
scope array of small imaging telescopes on La Palma in 
the Canary Isl ands run by an A rmenian-German-Spanish 
collaboration (Daum et al. 1997), and the Seven Telescope 
Array in Ut ah, which is opera ted by a group of Japanese 
institutions (Aiso et al. 1997). These techniques are rel- 
atively mature and the results from contemporaneous ob- 
servations o f the same source with different telescopes are 
consistent (Protheroe et al. 1997). Vigorous observing 
programs are now in place at all of these facilities. A vital 
observing threshold has been achieved whereby both galac- 
tic and extragalactic sources have been reliably detected. 
Many exciting results are anticipated as more of the sky 
is observed with this present generation of telescopes. 

The atmospheric Cherenkov imaging technique has now 
been adopted at a number of observatories whose proper- 
ties are summarized in Table pi. 

Above 30 TeV there are enough residual particles in the 
electromagnetic cascades that they can be detected at 
high mountain altitudes us ing arrays of particle detectors 
and fast wavefront timing ( png 1998 ). These arrays have 
large collection areas (>10,000m z ), good angular resolu- 
tion (^0.5°), moderate energy resolution (~ 100%), good 
duty cycle (100%) and large FOV (~2sr); however their 
ability to discriminate 7-rays from charged cosmic rays is 
severely limited. Despite the early promise of these ex- 
periments, which led to a considerable investment in their 
construction and operation, no verifiable detections have 
been reported by particle air shower arrays. 

3. GALACTIC SOURCES 

3.1. The Crab Nebula 

The Crab Nebula was first detected with a 37 pixel cam- 
era on the Whipple Observatory 10-m optical reflector in 
1989. This early and somewhat crude instrument yielded 
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Fig. 1. — The Whipple Observatory 10-m imaging atmospheric Cherenkov telescope. 



Table 1 
Operating ACIT Observatories c. 1999 May 



Group 


Countries 


Location 


Telescopc(s) 


Camera 


Threshold 


Epoch 








Number X Aperture 


Pixels 


(TeV) 


Beginning 


Whipple 


USA-UK-Irel. 


Arizona,USA 


10 m 


331 


250 


1984 


Crimea 


Ukraine 


Crimea 


6x2.4m 


6x37 


1 


1985 


SHALON 


Russia 


Tien Shen, Russia 


4m 


244 


1.0 


1994 


CANGAROO 


Japan-Aust. 


Woomera,Aust. 


3.8m 


256 


0.5 


1994 


HEGRA 


Gcrman-Armen.-Sp. 


La Palma, Sp. 


5x3m 


5x271 


0.5 


1994 


CAT 


France 


Pyrenees 


3 m 


600 


0.25 


1996 


Durham 


UK 


Narrabri,Aust. 


3x7m 


1x109 


0.25 


1996 


TACTIC 


India 


Mt. Abu, India 


10 m 


349 


0.3 


1997 


Seven TA 


Japan 


Utah.USA 


7x2 m 


7x256 


0.5 


1998 



a 9 a d etection with some 60 hours of integration on the 
source ( Wcckcs ct al. 198S ). The detection relied on dis- 
crimination of the 7-ray images from the much more nu- 
merous hadron background images. The possibility of a 
systematic effect with a new, and not yet proven, tech- 
nique could not be completely discounted (although nu- 
mero us tests were made for consistency). Nonetheless it at energies between 300 GcV and 50 TeV ( Hillas ct al 

requii ed the independent confirmation of the detection by 1998 Tanimori ct al. 1998b| ). To date there have been 



energy (100 MeV). 

Since then the Crab has been detected by eight inde- 
pendent groups using different versions of the atmospheric 
Cherenkov imaging technique (including one group in the 
Southern Hemisphere) . Some of these detections are listed 
in Table 0. The energy spectrum is now wel l determined 



other groups using different versions of the technique to 
really convince skeptics. The Whipple group subsequently 
detected the source at the 20 a level usin g an upgraded 
camera (109 pixels) ( Vacanti et al. 1991 ) and now rou- 
tinely detects the source at the 5-6 a level in an hour of ob- 
serva tion. The detected photon rate (about 2 per minute) 
is mo :e than that registered by EGRET at its optimum 



no positive detections reported by air shower array experi- 
ments using part icle detect ors (which operate at somewhat 
higher energies) ( Ong 1998 ) 



The simple Compton-synchrotron model (Gould 1965) 
has been updated to t ake account of a better under- 
standing of t he nebula ( |de Jager fc Harding 1992 ; Hillas 



et al. 1998); the measured flux is in good agreement 
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Table 2 
Flux from the Crab Nebula 



Group 



VHE Spectrum 
(10 -11 photons cm~ 2 s _1 ) 



Eth 
(TeV) 



(25(E/0.4TeV))-^ 4±u ^ 0X~ 

(3.2 ± o.7)(E/TcV)(- 2 ' 49±0 - 06 = t -" t±ao5 ^ t 0.3 

(2.7 ± 0.2 ± o.8)(E/TeV)- 2 - 61±ao6 ^ t±0 - 10 ^ t 0.5 

p. 7 ± 0.17 ± 0.40)(E/TeV)~ 2 - 57±0 - 14 ^ t±0 - 083 >- t 0.25 



Whipple (1991) a 

Whipple (1998) 6 

HEGRA (1999) c 

CAT (1998) rf 






acanti et al. 1991 



illas et al. 1998 
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Fig. 2 — The VHF, spe ctral energy distribution of the Crab Nebula compared with the predictions of a synchrotron self-Compton emission 
model (Hillas et al. 1998). 



with the predicted flux for a value of magnetic field 
(1.6±0.1 x 10~ 4 G) that is slightly lower than the equipar- 
tition value (Figure g) . Although this model is certainly 
simplistic given the structure now seen in optical images 
of the nebula, it shows that there is a viable mechanism 
that must work at some level. 

As in many other bands of the electromagnetic spec- 
trum, the Crab Nebula has become the standard candle 
for TeV 7-ray astronomy. Most importantly perhaps it is 
available as a steady source to test and calibrate the ACIT 
and can be seen from both hemispheres. Improvements in 
analysis techniques developed on Crab Nebula data have 
led directly to the detections of the AGNs discussed below. 

3.2. Supernova Remnants: Plerions 



The Crab Nebula is a somewhat unique object and hence 
one could not confidently predict what other supernova 
remnants might be detectable. The Crab is a member of 
that sub-class of supernova remnants known as plerions in 
which a bubble of relativistic particles is powered by a cen- 
tral pulsar. No other plerions have been seen by teles copes 
in the Northern Hemisphere ( Reynolds et al. 1993| ) (Ta- 
ble ||). However two have been detected in the Southern 
Hemisphere by the CANGAROO group. They first re- 
ported a detection o f PSR1706-44 at TeV energies in 1993 
( Kifunc et al. 1995 ) based on sixty hours of observation 
in the summer of 1992. PSR1706-44 is identified with a 
pulsar (of period 102 ms) and appears to be associated 
with a supernova remnant, possibly a plerion. At GeV 
energies it has a very flat spectrum. The energy spectrum 
is hard with a flux above 1 TeV of about 0.15 x 10~ n - 
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Table 3 
TeV Observations of Plerions 



Source Energy Flux/Upper Limit 
(GeV) (xlCT 11 ernes' 1 ) 



Group 



Crab Nebula 

PSR 1706-44 

Vela 
SS433 

3C58 
PSR0656+14 



400 

1000 
2500 
550 
550 
1000 



7.0 Whipple, ASGAT, HEGRA, TA 

Crimea, *Gamma, CANGAROO, CAT 

0.8 CANGAROO, Durham 

0.29 CANGAROO 

<1.8 Whipple 

<1.1 Whipple 

<3.4 Whipple 



There is no evidence that the signal is peri- 
odic. The detection has been confirmed by the Uni versity 
of Du rham grou p working in Narrabri, Australia ( |Chad 
wick 4t al. 1997|) . 



The CANGAROO group have also reported the detec- 
tion of a 6<r signal from the vicinity of the Vela pulsar 
(Yoshikoshi et al. 1997). The integral 7-ray flux above 



2.5 TeV is 2.5 x lO -1 ^ photons cm~ 2 s . Again there is 
no evidence for periodicity and the flux limit is about a 
factor of ten less than the steady flux. The signal is offset 
(by 0.14°) from the pulsar position which makes it more 
likely that the source is a synchrotron nebula. Since this 
offset position is coincident with the birthplace of the pul- 
sar it is suggested that the progenitor electrons are relics 
of the initial supernova explosion and they have survived 
because the magnetic field was weak. 

3.3. Pulsars 



The po wer spectra of most of the 7-ray pulsars (Thomp- 
son 1997) are extremely flat with maximum power often 
coming in the GeV energy range (see Figure ra). Because 
pulsar models often involve electrons with energies up to 
10 15 eV, it would come as no surprise if TeV 7-rays should 
emerge from the pulsar magnetosphere and be detected. 
Although there is, as yet, no established model for high 
energy 7-ray emission from pulsars, it appears that, in gen- 
eral, the detection of pulsed TeV 7- rays would favor outer 
gap (Romani 1996 ) over polar cap ( Daugherty fc Harding 
1982 ) models. This is because, in the latter models, the 
TeV 7-rays are attenuated by pair-production interactions 
with the intense magnetic fields near the pulsars. 

Sensitive upper li mits have been obtai ned for emis sion 
by the Crab pulsar (Lessard et al. 199E), Gcminga flAk-l 

ri ct| 



1993) and the Vela pulsar (Yoshikoshi 



erlof et al. 

1997[ ); in general these confirm the steepening of 
the spectra seen at 10 GeV energies. The radio pulsar 
PSR 1951+32 is particularly interesting because the power 
spectrum indicates that the maximum power occurs at en- 
ergies of at least a few GeV (Figure ||); in the EGRET 
measurements there is no evidence for a high energy cut- 
off. In fact, the flux continues to rise with energy up to 
the highest energy observation. Outer gap models suggest 
that the pulsar should be detectable at higher energies. 
Observations of PSR 1951+ 32 by the Whipple gro up re- 
ported only an upper limit flBrinivasan et al. 1997 ). This 
upper limit to the pulsed flux is two orders of magnitude 
below the flux extrapolated from the EGRET measure- 



ments. This represents the most dramatic turn-over in 
the spectrum of a 7-ray pulsar and hence puts the most 
severe constraints on the models. 



3.4. Supernova Remnants: Shell 

Supernova remnants (SNRs) are widely believed to be 
the sources of hadronic cosmic rays up to energies of ap- 
proximately Z x 10 14 eV, whe re Z is the nuclear charge of 
the particle (for a review see Jones et al. 1998J ). The ar- 
guments in support of this statement arc two-fold. First, 
supernova blast shocks are some of the few galactic sites 
capable of satisfying the energy required for the produc- 
tion of galactic cosmic rays, although even these must have 
a high efficiency, ~10% - 30% (e.g., Drury, Markiewicz & 
Volk 1989), for converting the kinetic energy of supernova 
explosions into high energy particles. Second, the model 
of diffusive shock accele ration 
1978| ; [Bell 19"7S| ; |Legage 



Blandford & Ostriker 



Cesarsky 1983), which provides 
a plausible mechanism for efficiently converting the explo- 
sion energy into accelerated particles, naturally produces 
a power-law spectrum of dN/dE ex E~ 21 . This is consis- 
tent with the inferred spectral index at the source for the 
observed local cosmic-ray spectrum of dN/dE ex E~ 2 - 7 , af- 
ter c orrecting for the effe cts of propagation in the galaxy 



(e.g., S wordy et al. 1990) 



The origins of cosmic rays cannot be studied directly 
because interstellar magnetic fields isotropize their direc- 
tions, except perhaps at the highest energies (> 10 18 eV). 
Thus, we must look for indirect signals of their presence 
from astrophysical sources. That SNRs accelerate elec- 
trons to high energies is well-established from observations 
of synchrotron emission in t he shells of SNRs at rad io and, 
more recently, X-ray (e.g., Koyama et al. 1995) wave- 
lengths. However, it is difficult to extrapolate from these 
data to inferences about the nature of the acceleration of 
hadrons in these same objects. 

Evidence of shock acceleration of hadronic cosmic rays 
in SNR shells could come from measurements of 7-ray 
emission in these objects. Collisions of cosmic-ray nuclei 
with the interstellar medium result in the production of 
neutral pions which subsequently decay into 7-rays. The 
7-ray spectrum would extend from below 10 MeV up to 
^1/10 of the maximum proton energy (;> 10 TeV), with 
a distinctive break in the spectrum near 100 MeV due to 
the resonance in the cross-section for 7r° production. As 7- 
ray production requires interaction of the hadronic cosmic 
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Fig. 3. — Power spectrum of 7-ray pulsars detected with EGRET (Thompson 1997) 



rays with target nuclei, this emission should be stronger 
for th ose SNR located near, or interacting with, dense tar- 



gets, mch as molecular clouds. The cosmic-ray density, 



and hence the associated 7-ray luminosity, will increase 
with time as the SNR passes through its free expansion 
phase, will peak when the SNR has swept up as much 
interstellar material as contained in the supernova ejecta 
(the Sedov phase) and gr adually decline thereafter (Drury, 
Aharonian & Volk 1994; |Naito fe Takahara 1994J ). Thus, 
7-ray bright SNRs should be "middle-aged." 

From the calculations of Drury et al. (1994), the lumi- 
nosity of 7-rays from secondary pion production may be 
detectable by the current generation of satellite-based and 
groun d-based 7-ray detectors, particularly if the objects 



identify the detected objec t with the SNR sh e ll. Because 
of this, embedded pulsars ( Brazier ct al. 199q ; |dc Jager fc| 
Mastichiadis 1997|; Harrus, Hughes fc Hclfand 1996) and 



an X-ray binary (|Kaaret et al. 199E) consistent with the 



positions of some of these EGRET sources have been sug- 
gested as alternative counterparts. In addition, significant 
background from the diffuse Galactic 7-ray emission com- 
plicates spectral measurements. To complicate matters 
further, with the detection of X-ray synchrotron radia- 
tion from SNR shells, the possibility that 7-rays could be 
produced via inverse Com pton scattering of ambient soft 



are lo :ated in a region of relatively high density in the in- 



photons has been r ealized (Mastichiadis & de Jager 1996 
Mastichiadis 1996). Bremsstrahlung radiation may also 
be a significant source of 7- rays at MeV-GeV energies ( pq 
Jager & Mastichiadis 1997; Gaisser, Protheroe & Stanev 



terstellar medium. In support of this hypothesis, EGRET 
has detected signals from several regions of th e sky con- 
sisten t with t h e positions of shell-typ e SNRs ( Bturncr fc 
De rnier 1995; [Esposito ct al. 1996 ; Lamb & Macomt 
19971 jjaffe etal. 19971). However, the EGRET detections 



1998) 



3.4.1. 



alone are not sufficient to claim the presence of high energy 
hadronic cosmic rays. For instance, the relatively poor an- 
gular resolution of EGRET makes difficult to definitively 



VHE j-ray observations 

Measurements of 7-rays at very high energies may help 
resolve the puzzle of the 7-ray emission from the EGRET- 
detected sources. VHE 7-ray telescopes have much bet- 
ter angular resolution than EGRET, reducing the source 
confusion associated with any detection. Also, because 
the diffuse Galactic 7-ray emission has a relatively steep 
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Table 4 
Observations of shell-type supernova remnants 



Object 
Name 


Observation 
Time 
(min.) 


Energy 
(TeV) 


(io- 


Integral 
Flux a 

-1 1 —2 
li cm z 


B^) 


Rcf. 


Tycho 


867.2 

1076.7 

678.0 

360.1 

468.0 

560.0 

2820.0 

140.0 

2040.0 


>0.3 
>0.3 
>0.5 
>0.3 
<0.3 
>0.3 
>0.5 
>0.3 
>1.7 


0.46 ± 0.6 stat i 


<0.8 
<2.1 

<1.9 C 
<3.0 
<3.6 
<2.2 

<l.l c 
<6.4 


Buckley et al. 1998 




IC 443 6 


Buckley et al. 1998 






Hess et al. 1997 




W44 b 


Buckley et al. 1998 




W51 


Buckley et al. 1998 




7-Cygni 6 
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a Upper limits from Buckley et al. (1998) are at the 99.9% confidence level and those from 

Hess et al. (1997) are at the 3 a confidence level. 

6 Associated with an EGRET source. 

c Upper limits were converted from fractions of >500GeV Crab flux using the measured 

Crab flux of Hillas et al. (1998). 



27 ( |Hunter et al. 1997| ), com- 
~ E - ^ 1 spectrum of 7-rays from 



spectrum, oc E — E 

pared with the expected 

secondary pion decay, contamination from background 7- 

ray emission should be less in the VHE range. Thus, in 

recent years, searches for emission from shell-type SNRs 

have been a central part of the observation program of 

VHE telescopes. 

The Whipple Collaboration has published the results 
of observations of six shell- type SNRs (IC443, 7-Cygni, 
W44, W51, W63, and Tycho) selected as strong 7-ray 
candidates based on their radio properties, distance, small 
angula r size, and possible association with a molecular 
cloud ( |Buckley et al. 1998| ; [Hess ct al. 1997| ). The small 
angular size was made a requirement due to the limited 
field of view (3° diameter) of the Whipple telescope at 
that time. VHE telescopes can also detect fainter 7-ray 
sources if they are more compact, because they can reject 
more of the cosmic-ray background. IC443, 7-C ygni, and 
W 44 are a lso associated with EGRET sources ( Esposito] 



et al. 1996). Despite long observations, no significant ex- 



cesses were observed, and stringent limits were derived on 
the VHE flux (see Table @). 

In contrast to the upper limits derived by the north- 
ern hemisphere telescopes, the CANGAROO Collabora- 
tion has recently reported evidence for TeV 7-ray emis 



sion from the shell- type SNR, SN 1006 (Tanimori ct al 
1998a). Observations taken in 1996 and 1997 indicate a 
statistically significant excess from the northeast rim of 
the SNR shell (see Figure 0) . The position of the excess 
is consistent with the location of non-thermal X rays de- 
tected by the ASCA experiment (Koyama ct al. 1995). If 
this object is confirmed as a TeV 7-ray source, it repre- 
sents the first direct evidence of acceleration of particles 
to TeV energies in the shocks of SNRs. 

3.4.2. Implications of the j-ray observations 

If the EGRET detections do indicate the presence of 
7-rays produced by secondary pion decay, the measured 



flux can be compared with the VHE upper limits. These 
VHE upper limits and the EGRET measurements are com- 
pared to the predicted fl uxes from the model of Drury et 
al. (1994) in Figure | flBucklcy ct al. 1998| ). The solid 
curves are normalized to the integral >100McV flux de- 
tected by EGRET assuming a source cosmic-ray spectrum 
of E -2 - 1 , so they assume that the EGRET emission is en- 
tirely due to cosmic-ray interactions and that the cosmic- 
ray spectrum at the source has the canonical spectrum. 
In the cases of 7-Cygni, IC443, and W44, the Whipple 
upper limits lie a factor of ~25, 10, and 10, respectively, 
below the model extrapolations and require either a spec- 
tral break or a differential source spec trum steeper than 
E -2.5 for 7 _Cygni and E~ 2 - 4 for IC 443 ( [Buckley 199^ ). In 
addition, the measured EGRET spectra, while consistent 
with a spectral index of about 2, do not show the flat- 
tening of the 7-ray spectrum near 100 MeV which would 
be expected if the 7-rays result from secondary pion de- 
cay. Gaisser et al. (1998) performed multiwavelength fits 
to the EGRET and Whipple results and concluded that 
if the EGRET detections are truly from the shells of the 
SNRs, the EGRET data must be dominated at low ener- 
gies by electron bremsstrahlung radiation, but the source 
spectrum at high energies must still be relatively steep 
( ~E~ 2,4 ) to accoun t for the Whipple upper limits (cf., 



Buckley et al. 1998) 



If, on the other hand, the > 100 MeV emission results 
from some other emission mechanism, such as the inter- 
actions of high energy electrons accelerated by embedded 
pulsars or in the SNR shocks, the EGRET results should 
not be compared to the Whipple data. Instead, the Whip- 
ple data must be considered alone in the context of the 
secondary pion decay models. This possibility was also 
investigated by Buckley et al. (1998) and is shown by the 
dashed curves in Figure ||. The dashed curves represent a 
conservative estimate of the range of allowable parameter 
values for the Drury et al. (1994) model, without reference 
to the EGRET detections. For this comparison, there is 
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Fig. 4. — Left: contour map of the statistical significance of the excess emission from SN 1006 as observed with the CANGAROO telescope 
in 1996. Right: the same plot for data taken in 1997. The solid lines indicate the contour map of non-thermal X-ray emission detected with 
ASCA. The dashed circles indicate the angular resolution of the CANGAROO telescope. Figure from Tanimori et al. (1998a). 



still room for the models to work in these objects. How- 
ever, the upper limits in some (e.g., IC443) are beginning 
to strain the limits of the available parameter space. It will 
take more sensitive measurements with future telescopes 
to fully span the allowable parameter space and see if we 
need to reconsider our assumptions about the sources of 



cosmi r, lays witliiii urn ualaAV. 

Th e TcV cmiaaion from SN 1006 detected by the CAN 
GAROO group also does not require the presence of 
hadronic cosmic rays. In fact, the most common expla- 
nation for the detected emission is inverse Compton scat- 
tering _of_cl_C£trOTis^thco^^ pho- 
tons ( [Reynolds 1996] ; [Mastichiadis k de Jager 1996J ). The 



its fo 



that Lin 



mam Higmiicnis lot tins aiu mat inu emission is iciiieieu 
on on e of the rcgiona where the synchrotron cmiaaion waa 
detected with ASCA and the lack of evidence for a nearby 
molecular cloud needed to boost the TeV 7-ray flux to a de- 
tectable level. Under the assumption that the emission is 
from inverse Compton scattering, Tanimori et al. (1998a) 
have combined their data with the ASCA results to derive 
an estimate of 6.5 ± 2 fiG for the magnetic field within 
the SNR shell. The observations also provide an upper 
limit on the acceleration time. Thus, the TeV observations 
provide previously unknown parameters for models of the 
shock acceleration in SNRs. Unfortunately, the possibility 
of inverse Compton emission from the shells of SNRs also 
confuses the issue for using 7-rays as a probe of cosmic-ray 
acceleration in SNRs. Future measurements will need to 
provide accurate spectra and spatial mapping of the 7-ray 
emission from SNRs in order for the source of the 7-ray 
emission to be unambiguously resolved. 

4. EXTRAGALACTIC SOURCES 



4.1. Active Galactic Nuclei 

In recent years, high energy 7-rays have come to play an 
important role in the study of AGNs. Before the launch 
of CGRO in 1991, the only known extragalactic source of 
high energy 7-rays was 3C 273 which had been detected 



with the COS-B satellite 20 years ago (Swanenburg et al 



1978). The EGRET detector on the CGRO has identified 



more tha n 65 AGNs which em it 7-rays at energies above 
100 McV ( Hartman et al. 1999| ), and a substantial fraction 
of those sources which remain unidentified in the EGRET 
catalog are likely to be AGNs as well. In addition, the 
Whipple Observatory 7-ray telescope has disc overed three 



AGN s which emit at ene r gies above 300 GeV ( Punch et al 
1992; Quinn et al. 1996; Catanese et al. 1998]) and there 



are recent detections of two other AGNs with Cherenkov 



telescopes (Chadwick et al. 1999; Neshpor et al. 1998) 



During flaring episodes, the 7-ray emission can greatly ex- 
ceed the energy output of the AGNs at all other wave- 
lengths. Thus, any attempt to understand the physics of 
these objects must include consideration of the 7-ray emis- 
sion. 

All of the AGNs detected in high energy 7-rays are 
radio-loud sources with the radio emission arising primar- 
ily from a core region rather than from lobes. These types 
of AGNs are often collectively referred to as "blazars," 
and include BL Lacertae (BL Lac) objects, flat spectrum 
radio- loud quasars (FSRQs), optically violent variables, 
and superluminal sources. The emission characteristics of 
blazars include high polarization at radio and optical wave- 
lengths, rapid variability at all wavelengths, and predom- 
inantly non-thermal emission at most wavelengths. The 
emission from blazars is believed to arise from relativistic 
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Fig. 5. — Whipple Observatory upper limits (W) shown along with EGRET integral fluxes (E) and integral spectra. Also shown are 
CASA-MIA upper limits (CM) from Borione et al. (1995), Cygnus upper limits (C) from Allen et al. (1995), and the AIROBICC upper 
limit from Prosch et al. (1996). The solid curves indicate extrapolations from the EGRET integral data points at 100 MeV (indicated by the 
triangles). The dashed curves are estimates of the allowable range of fluxes from the model of Drury et al. (1994). Figure from Buckley et 
al. (1998). 



jets oriented at small angles to our line of sight. If so, 
the observed rad iation will be strongly a mplified by rela- 
tivistic beaming (Blandford & Rees 1978). Direct evidence 
for relativistic beaming of the radio emission comes from 
Very Long Baseline Interferometer (VLBI) observatio ns of 
appar ent super luminal motion in many blazars (e.g., |Vcr- 
meulen & Cohen 1994). The rapid variability and high 



luminosities of the detected 7 -ray sources imply that the 
7-rays are also beamed (see § 4.1.3) ) . 

There is a growing consensus that blazars are all the 
same type of o bject, perhaps differ i ng only in intrinsic lu- 
minosity (e.g., Fossati et al. 1998; phisellini et al. 1998) 



or some combination of luminosity and viewing angle (e.g., 
Gcorganopoulos & Marscher 1998). However, for this work 



we wi ll continue the practice of referring to BL Lac objects 



and FSRQs as distinct objects: BL Lac objects are those 
blazars which have optical emission lines with equivalent 



width < 5 A and FSRQs are the remaining blazars. As 
we will see, this distinction may be important in explain- 
ing why only BL Lac objects are detected at very high 
energies. 

The spectral energy distribution of blazars appears to 
consist of two parts. First, a low energy component ex- 
hibits a power per decade distribution that rises smoothly 
from radio wavelengths up to a broad peak in the range 
spanning infrared (IR) to X-ray wavelengths, depending 
on the specific blazar type, above which the power output 
rapidly drops off. Second, a distinct, high energy compo- 
nent, which does not extend smoothly from the low energy 
component, is often seen. It typically becomes apparent 
in the X-ray range and has a peak power out put in the 
7-ray range betw een ~lMeV and ITeV (e.g., von Mon 



tigny et al. 1995), again depending on the specific blazar 
type. When plotted as E 2 dN/dE (or equivalently v¥ v ) the 
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determined by EGRET (Thompson et al. 19951). The position 
of Mrk421 is indicated by the cross, figure from Buckley et al. 
(1996). 
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Fig. 7. — Two-dimensional plot of the VHE 7-ray emission from 
the region around Mrk501. The gray scale is proportional to the 
number of excess 7-rays and the solid contours correspond to 1 a 
levels. The position of Mrk 501 is indicated by the cross. 



spectral energy distribution shows a two-humped shape, 
thou gh some objects sh ow evidence of a third component 



(e.g., Kubo et al. 1998) 



Though there is no general consensus on the origin of 
these emission components, it is generally agreed that 
the low energy component arises from incoherent syn- 
chrot ron emission by relativist ic electrons within the jet 



(e.g., Blandford fc Rees 1978|). This is supported most 



strongly by the high-level, variable polarization observed 
in these objects at radio and optical wavelengths. The 
origin of the high energy emission is a matter of great 
interest. There are many variations of the models and 
here we only briefly mention a few which are most often 
invoked to explain the 7-ray emission. The most pop- 
ular models at this time are those in which the 7-rays 
are produced through inverse Compton scattering of low 
energy photons by the same electrons which produce the 
synchrotron emission at lower e nergies. Syn c hrotron sclf- 



Compto n (SSC) emission (e.g., Konigl 1981 ; Maraschi et 
al. 19B3; Bloom & Marschcr 1996), in which the seed pho- 
tons for the scattering are the synchrotron photons already 
present in the jet, must occur at some level in all blazars, 
but models in which the 7-ray emission arises predomi- 
nantly from inverse Compton scattering of seed photons 
which arise outside of the jet, either directly from an ac- 
cretion disk (Dermer, Schlickeiser & Mastichiadis 1992) or 
after being re-processed in the broad-line region or scatter- 
ing off thermal plasma (Sikora, Begelman & Rees 1994), 
appear to fit the observations satisfactorily as well. An- 
other set of models proposes that t he 7-rays are prod uced 



by proton-initiat ed ca scades (e.g., Mannheim 1993). As 



we will show in §|4.1.3, the 7-ray observations strain both 



types of models, but do not, at present, rule any out. 



In the remainder of this section we discuss the status of 
VHE observations and t he emission characteristics of the 
detected objects (§ |4.1.l| ), the results of m ulti-wavelength 
campaigns on the detected objects (§4.1.2) which are the 
best probe of the physics of blazars, and finally briefly 
discuss some of the implications of these observations on 
our understanding of the physics of the blazars and on the 
models which purport to explain them. 

4.1.1. Observational Status and Emission Characteristics 

The BL Lac object Markarian421 (Mrk 421, z = 0.031) 
was detected as the first extragalactic source of VHE 7- 
rays in 1992 using the W hipple Observatory 7-ray tele- 
scope ( Punch et al. 1992| ). A two-dimensional image of 
the emission from Mrk 421 is shown in Figure pi Although 
Mrk 421 had previously been par t of an active program of 
observing extragalactic sources ( pawley et al. 1985 ) by 
the Whipple Collaboration, the observations which led to 
the detection of Mrk 421 at TeV energies were initiated in 
response to the detection of several AGN by the EGRET 
experiment. The initial detection indicated a 6 c excess 
and the flux above 500 GeV was approximately 30% of the 
flux of the Crab Nebula at those energies. Mrk 421 has 
been confirmed as a source of VHE 7-rays by the HEGRA 
Collabor ation (Petry et a l. 1996 ), the Telescope Array 
Project ([Aiso et al. 1997|), andthe SHALON telesc ope 



( Sinitsyna et al. 1997] ) and, as discussed in §4.1.2 be- 
low, multi-wavelength correlations have confirmed that the 
VHE source is indeed Mrk 421 and not some other object. 
With the successful detection of Mrk 421, the Whipple 
Collaboration initiated a search for VHE emission from 
several other blazar-type AGNs, concentrating at first on 
those objects detected by EGRET, but also spending a 
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Fig. 8. — Daily VHE 7-ray count rates for Mrk421 during 1995. Modified Julian Day 49720 corresponds to 1995 January 3. Figure from 
Buckley et al. (1996). 



substantial amount of time observing radio-loud blazars 
which were not detected by EGRET. This broad ap- 



proac 1 led to the detection of the BL Lac object Mrk 501 
(z = 0.034) by the Whipple Collaboration in 1995 flQuinrj 



300 GeV at the la level by combining o bservations from 
1996 and 1997 by the Durham group (|C had wick et al 



1999| ). The flux was approximately 40% of the VHE flux 



et al. 1996). A two-dimensional mar 



of the VHE emis- 
Because Mrk 501 



sion from Mrk 501 is shown in Figure u\ 
had not been detected as a significant source of 7-rays by 
EGRET, this was the first object to be discovered as a 7- 
ray source from the ground. Hence, VHE 7-ray astronomy 
was e stablished as a legitimate channel of astronomical in- 



of the Crab Nebula and corresponded to an active X-ray 
emission period. PKS 2155-304 is an EGRET source with 
an average flux at E > 100 MeV comparable to that of 
Mrk 421. Finally, the BL Lac object 3C66A (z = 0.444) 
has been reported as a source of >900GeV 7-rays based 
on a 5 a excess seen in observations in 1996 by 7-ra y tele- 



scopes at the Cr imean Astrophysical Observatory ( Nesh 



vestig itions in its own right, not just an adjunct of high por et al. 1995). The average flux during these obser 



energy observations from space. The flux of Mrk 501 dur 
ing 1995 was, on average, 10% of the VHE flux of the 
Crab Nebula, making it the weakest detected source of 
VHE 7-rays. Mrk 501 was confirmed as a source of VHE 



7-rays in 1996 by the HE GRA telescop es ( Bradbury et al 
19970 the C AT telescope ([Punch 1997|) , the Telescop e Ar- 



ray Project flHayashida et al. 1998| )7and TACTIC (ghat 
19970 

In addition to the confirmed detections of Mrk 421 and 
Mrk 501, three other objects have recently been reported 
as sources of VHE 7-rays, but remain to be verified by de- 
tections from independent 7-ray telescopes. The BL Lac 
object 1ES 2344+514 (z = 0.044) was detected at ener- 
gies about 350 GeV by the Whipple Observatory in 1995 
( Catanese ct al. 1998f ). Most of the emission comes from 
a single night, December 20, in which a flux of approxi- 
mately half that of the Crab Nebula was detected with a 
significance of 6 a. Other observations during that year 
revealed an excess of 4 er; subsequent observations have 
yielded no significa nt signal. 1ES 2 344+514 is not de- 
tected by EGRET ( Thompson 1996 ) so if this detection 
is confirmed it is another instance of a 7-ray source be- 
ing first detected by a ground-based telescope. PKS 2155- 
304 (z = 0.117), often considered the archetypical X-ray 
selected BL Lac object, was detected at energies above 



vations was approximately 120% of the flux of the Crab 
Nebula at these ener gies. 3C 66A is an EGRET source 
( |Hartman et al. 1999| ). 

Extreme variability on time-scales from minutes to years 
is the most distinctive feature of the VHE emission from 
these BL Lac objects. Variability in the emission is a sur- 
prising feature in some respects because it implies a small 
emission region. If low energy photons (e.g., infrared, op- 
tical, and ultraviolet) are produced in the same region, 
the VHE photons would pair produce with these photons 
and would not escape. Also, if the variability occurs near 
the base of the jet, there is likely to be considerable ambi- 
ent radiation present which can attenuate the 7-ray signal. 
This opacity problem is reduced considerably if the emis- 



sion is beamed towa rd us (e.g., Dermcr fc Gehrcls 1995 



Buckley et al. 1996| ), and this has been one of th e mai n 



arguments for 7-ray beaming in these objects (see § 4.1.3). 

The first clear detection of flaring activity in the VHE 

emission of an AGN came in 1994 observations of Mrk 421 



by the Whipple Collaboration ( Kerrick et al. 1995a ) where 
a 10-fold increase in the flux, from an average level that 
year of approximately 15% of the Crab flux to approx- 
imately 150% of the Crab flux, was observed. Subse- 
quent analysis of Mrk 421 data indicated evidence for less 
prominent episodes of variability during 1992 and 1993 as 
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Fig. 9. — Lightcurves of two flares observed from Mrk421 by the Whipple Collaboration on 1996 May 7 (a) and May 15 (b). The time axes 
are shown in coordinated universal time (UTC) in hours. For the May 7 flare, each point is a 9-minute integration; for the May 15 flare, the 
integration time is 4.5 minutes. Figure from Gaidos et al. (1996). 



well ( Bchubncll ct al. 1996 ). This suggested that vari- 
ability could be present on a fairly frequent basis in the 
VHE emission. In order to characterize this variability, 
the Whipple Collaboration began systematic monitoring 
of Mrk421 in 1995 which continues to the present day. 
The observations of Mrk421 in 1995, shown in Figure g, 
revealed several distinct episodes of flaring activity, like 
in previous observations, but, perhaps more importantly, 
indicated that the VHE emission from Mrk421 was best 
characterized by a succession of day-scale or shorter flares 
with a baseline emission level below the sens itivity limit 
of the Whipple detector ( Buckley ct al. 1996J) . The time- 
scale of the flaring is derived from the fact that, for the 
most part, the flux levels measured each night varied fairly 
randomly with no evidence of a smooth pattern. Thus, 
thoug h no significant intra-night variability was discerned 



in these observations, it seemed clear that it could occur. 



The hypothesis that the VHE emission of Mrk421 could 
flare on sub-day time-scales was borne out in spectacular 
fashion in 1996, with the observat ions of two flares by the 
Whipple Collaboration (Figure O; Gaidos et al. 1996). In 



the Crab, at which point observations had to stop because 
of moonrisc. This flux is the highest observed from any 
VHE source to date. The doubling time of the flare was 
~1 hour. Follow-up observations on May 8 showed that 
the flux had dropped to a flux level of «30% of the Crab 
Nebula flux, implying a decay time-scale of <1 day. The 
second flare, observed on May 15, although weaker was re- 
markable for its very short duration: the entire flare lasted 
approximately 30 minutes with a doubling and decay time 
of less than 15 minutes. These two flares are the fastest 
time-scale variability, by far, seen from any blazar at any 
ry-ray energy. 

Systematic observations of Mrk501 sensitive to day- 
scale flares have been conducted si nce 1995 with the W hip- 
ple Observatory 7-ray telescope ( Quinn ct al. 1999|) and 
since 1997 wi th the te lescopes of th e HEGRA ( Aharonian 
et al. 19991), CAT (|Punch 19971) , and Telescope Array 
flHayashida et al. 199S ) collaborations. The results of 



the first flare, observed on May 7, the flux increased mono- 
tonically during the course of ~2 hours of observations, be- 
ginning at a rate twice as high as any previously observed 
flare and reaching a counting rate wlO times the rate from 



these observations indicate a wide range of emission lev- 
els (Figure |l0[) and some very interesting similarities and 
differences with the VHE emission from Mrk421. The ob- 
servations in 1995 indicate a flux which is constant, with 
the exception of one night, MJD 49920, when the flux was 
approximately 4.6 a above the average (approxim ately 5 
times the flux during the remainder of the season) (Quinn 
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Fig. 10. — VHE 7-ray lightcurve for Mrk501 as observed with the Whipple telescope between 1995 and 1998 at energies above 350 GeV. 
Fluxes are expressed as fractions of the Crab Nebula flux above 350 GeV. Figure adapted from Quinn et al. (1999). 



et al. 1996; Quinn ct al. 1999J ). Observations in 1996 
by the Whipple Observatory show that the average flux of 
Mrk501 had increased to approximately 20% of the Crab 
Nebula flux above 300 GeV indicating a two-fol d increase 
in th e ave rage flux over the 1995 observations (Quinn ct 
al. l |9 9 9| ) . HEGRA observations indicated an average 
flux of approximately 30% of the Crab flux above 1.5 TeV, 
perhaps indicating a harder emission spectrum than that 
of the Crab Nebula. The Whipple observations show no 
clear flaring episodes but the probability that the aver- 



age monthly flux levels are drawn from a distribution with 



(whereas previous observations had never revealed a flux 
>50% of the Crab flux). Also, the amount of day-scale 
flaring increased and, for the first time, significant hour- 
scale variations were seen. Two clear episodes of hour-scale 
variability were detected with the Whipple Observatory 
telescope (Figure |ll|) and a search for intraday variabil- 
ity revealed several other nights which, considered alone, 
would not have been considered significant but, when com- 
bined, indicated frequent intra-day variability whi ch was 



just below t he sensitivity of the Whipple telescope ( Quinn 
et al. 199S|). Analysis of data from the HEGRA (lAharo- 



a con stant flux level is 3.6x10 b , clearly indicating that nian et al. 1999a ) and Telescope Array ( Hayashida et al 



al. 1!'9S). There is no significant evidence for day-scale 
variations within each month in 1996. 

In 1997, the VHE emission from Mrk 501 changed dra- 
matically. After being the weakest known source in the 
VHE sky in 1995 and 1996, it became the brightest, with 
an average flux greater than that of the Crab Nebula 



the er iission is varying on at least month-scales (|Quinn et] 1998 ) projects revealed no statistically significant intra- 



night variations, but the two nights in the HEGRA data 
with the smallest statistical probability of having constant 
emission are the same nights seen to have significant vari- 
ability in the Whipple data. 

Perhaps the most important aspect of the observations 
of Mrk 501 in 1997 was that, for the first time, Chcrcnkov 
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Fig. 11. — Very high energy 7-ray light curves of Mrk501 for the two nights in 1997 which show significant intra-night variability. Figure 
from Quinn ct al. (1999). 



telescopes other than the Whipple telescope consistently 
detected a significant excess from Mrk501 on a nightly 
time-scale. This permitted more complete VHE light- 
curves to be obtained (which it is expected will eventu- 
ally lead to a better understanding of the VHE emission 
from blazars) and also provided confirmation that differ- 
ent VHE telescopes could obtain consistent results from a 
variable source (see Figure ^2|) . 

In addition to the establishment of the flaring itself, 
the Telescope Array Collaboration performed a periodic- 
i ty search with their VH E observations of Mrk501 in 1997 
( Hayashida et al. 1998 ). They show evidence for a quasi- 
periodic signal in the data which has a period of approx- 
imately 12.7 days. This is disturbingly close to half the 
lunar cycle. If real this would be an extraordinary result, 
given the very short time-scale of the quasi-periodicity. 

The only published results on observations of Mrk 501 in 



1998 are those of the Whipple Observatory ( Quinn ct al. 
19990 The average emission level was approximately 30% 
of the Crab flux but there was considerably more variabil- 
ity in the emission than in 1996 or 1995. Two distinct, 
very high flux flares were observed, one with the highest 
flux (approximately 5 times the Crab flux) ever observed 
from Mrk 501 with the Whipple telescope. The monthly 
average flux was also variable, with three months show- 
ing emission levels similar to the 1995 flux, approximately 



10% of the Crab Nebula. 

A natural question to ask about the variability in 
Mrk 501 is whether the degree of variability seen changes as 
a function of the mean flux level. That is, whether Mrk 501 
is really more variable when its average flux is higher or 
whether it is an artifact of the telescopes being more sen- 
sitive to variations when the average flux is higher? To 
test this, Quinn et al. (1999) performed simulations to 
see if the day-scale variability observed with the Whip- 
ple telescope in 1997 (when the average VHE flux was 1.3 
times that of the Crab) would have been detectable in 1996 
and 1995 (when the average VHE flux was 20% and 10% 
that of the Crab, respectively) and also tested whether 
the month-scale variations in 1997 and 1996 would have 
been detectable in 1995. Their simulations indicate that 
the day-scale flaring in 1997 would have been detectable in 
the 1996 data, but not the 1995 data, and the month-scale 
variations in 1996 would have been detectable in 1995, 
while the 1997 month-scale variations would not have been 
detectable. Thus, it appears that the higher state emis- 
sion levels have different variability characteristics than 
the lower emission levels. 

The other unconfirmed sources are, if they are in- 
deed sources, also variable emitters of VHE 7-rays. 
1ES 2344+514 was only detected with high statistical sig- 
nificance on one night, but has never been detected since, 
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Fig. 12. — Very high energy 7-ray observations of Mrk501 with the Whipple (filled circles, E>350GeV), CAT (open squares, E>250GeV), 
HEGRA (filled stars, E>lTeV), and Telescope Array (open triangles, E>600 GeV) telescopes between 1997 February and October. Fluxes 
are shown in units of the Crab Nebula flux for the energy threshold of each telescope. The numbers on the horizontal axes for each plot 
indicate the Modified Julian Days during that month. Data are from Quinn et al. (1999), Aharonian et al. (1999a), Djannati-Atai et al. 
(1999), and Hayashida et al. (1998). 



with lux limits of approximately 10% of the Crab flux et al. 1997; Krennrich et al. 1999). Analysis of the spec- 



( ICatancsc et al. 1998| [Aharonian et al 1999b|). PKS2155 
304 has also been claimed to be variable (Chadwick et al 



tra obtained from observations of flares on 1996 May 7 
and 15 and observations of high state emission taken at 
large zenith angles in 1995 June indicate that, within the 
statistical uncertainties, the spectra are all cons istent with 
a simple pow er law spectrum: dN/dE oc E -25 (Krennrich 
tory telescope flKcrrick et al." 19951) and the HEGRA et al. 1999). When combined, these three data sets are 



19991)] although the statistical probability that the emis- 
sion is constant is not quoted. Finally, 3C 66A must be 
variable, or else observations with the Whipple Observa- 



telescope ( |Aha.rouian et al. 1999b ) would have easily de- 
tected this object at the flux level quoted by the Crimean 
group. 

The high flux VHE emission from Mrk421 and Mrk501 
has permitted detailed spectra to be extracted. Accurate 
measurements of the VHE spectrum are important for a 
variety of reasons. First, the shape of the high energy spec- 
trum is a key input parameter of AGN emission models, 
particularly as it relates to the MeV-GeV measurements 
by EGRET. Second, how the spectrum varies with flux, 
compared to longer wavelength observations provide fur- 
ther emission model tests. Third, spectral features, such 
as breaks or cut offs, can indicate changes in the primary 
particle distribution or absorption of the 7-rays via pair- 
production with low ene rgy photons at the source or in 
intergalactic space (see § 4.2). 

For Mrk421, the only detailed spectra published at this 
time come from observations of high state emis sion with 
the Whipple Observatory telescope (Figure O; Zwccrink 



consistent w ith a simple power law spectrum for Mrk421 
of the form (Krennrich et al. 1999J) 



^| (250GeV-10TcV) K E - 2 - 54±003 — ±010 — 

where E is in TeV. 

Observations of Mrk421 in 1997 and 1998 with the 
HEGRA system of Cheren kov telescopes reveal a signif- 
icantly different spectrum ( |Aharonian et al. 1999c ) 



dN 
dE 



(500 GeV - 7 TeV) oc E- 309±007 =^ t±0 



10sy 



than observed with the Whipple telescope. The emission 
level for the HEGRA observations was approximately 0.5 
times the Crab flux, much lower than the fluxes (1 - 10 
times the Crab flux) used in the Whipple observations. 
This may indicate that the spectrum in Mrk421 becomes 
softer with decreasing flux. However, HEGRA observa- 
tions show no evidence of variability between observations 
at fluxes above the Crab flux and those between one-sixth 
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Fig. 13. — VHE 7-ray spectra of Mrk421 (open stars) and 
Mrk 501 (filled circles! as measured wi th the Whipple Observatory 
telescope (from Krennrich et al. 1999). The solid line through the 
Mrk 421 points indicates the best-ht spectrum to these data and 
the dashed line through the Mrk 501 points is the best-fit spectrum 
for those data. 
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the HEGRA array (filled circles) and the Whipple telescope (open 
circles). The dashed line indicates the power law plus exponential 
cut-off spectrum fit to the HEGRA data. Figure from Konopelko 
(1999). 



and one-half the Crab flux and the Whipple results show 
no variations in spectral index despite using observations 
spanning a 10-fold range of fluxes. Further studies may 
help resolve these differences. 

As with the studies of the variability of its VHE flux, the 
high state emission detected from Mrk 501 in 1997 allowed 
detailed spectra to be derived by several experiments, per- 
mitting studies of the time-dependence of the spectra and 
providing all-important cross-checks of the methods used 
to derive energy spectra. Because of the rapid variabil- 
ity of the emission, again the normalization of the spectra 
are not of fundamental importance, except perhaps in the 
c ontex t of multi- wavelength studies which are discussed in 
§ EO below. 

The most detailed energy spectra published at this time 
come fr om Whipple observat i ons between 250 GeV and 
12 TeV ( |Samuelson et al. 199S| ; [Krennrich et al. 1999|) and 
HEG RA d ata spanning 500 GeV to 20 TeV flAharonian et 
1999q). The Telescope Array Collaboration has also 



and the HEGRA spectrum is: 



al 



derived a spectrum o ver a slightly narrower energy range 
(600 GeV to 6.5 TeV) flHayashida est al. 1998Q . A search for 
variability in the spectrum re vealed no significant changes 
in spe ctrum w ith flux or time (]5amuelson 199S ; Aharoniar] 



et al. 1999a), allowing large data sets to be combined to 



derive very detailed energy spectra spanning large ranges 
in energy. The spectra derived by Whipple and HEGRA 
deviate significantly from a simple power law. For Whip- 
ple, the \ 2 probability that a power law is consistent with 
the measured spectrum is 2.5 x 10~ 7 . This is the first 
significant deviation from a power law seen in any VHE 7- 
ray source and any blazar at energies above lOMeV. The 
Whipple spectrum is: 



— ex E _2 - 22±0 - 045tat±0 - 05syBt_(0 ' 47±0 - 07stat ' O910 (' E ^ 
dE 



dN 
dE 



x 



g-1.92±0.03 8t at±0.20 sys t 



* exp 



E 



6.2±0.4 stat (-1.5 +2.9) 



syst 



where E is in units of TeV. The form of the curvature term 
in the spectra has no physical significance as the energy 
resolution of the experiments is not sufficient to resolve 
particular spectral models. The Whipple spectral form 
is simply a polynomial expansion in logE v. log(dN/dE) 
space. The HEGRA form was chosen presumably because 
attenuation of the VHE 7-rays by pair-production with 
background IR photons could produce an exponential cut- 
off. In fact, the Whipple and HEGRA data are completely 
consistent with each other as shown in Figure [lj. The 
Telescope Array Collaboration derived a spectrum which 
is well fit by a simple power law (dN/dE oc E~ 2 - 5±0 - 1 ). 
The data from this spectrum are also consistent with the 
Whipple and HEGRA spectra. 

4.1.2. Multi-wavelength observations 

Some of the most exciting results on VHE 7-ray sources 
have come through observations where several telescopes 
operating at different wavelengths simultaneously moni- 
tor the activity in a blazar. These multi-wavelength cam- 
paigns have involved the larger astronomical community in 
the study of VHE sources and served the subsidiary pur- 
pose of confirming the source identifications of the VHE 
7-ray emitting blazars. 

The first evidence of correlated variability between VHE 
7-rays and lower energy emission came from a multi- 



wavelength campaign on Mrk 421 in 1 994 April/May (Ma- 
comb et al. 1995, Macomb et al. 1996| ). Observations were 



conducted with t 



qc Whipple telescope, EGRET, the Ad- 



vanced Satellite for Cosmology and Astrophysics (ASCA) 
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Fig. 15. — Lightcurves for observations of Mkn421 in 1994 May by ASCA (top), Whipple (middle), and EGRET (bottom). Figure from 
Takahashi, Madejski, and Kubo (1999). 



in X rays, the International Ultraviolet Explorer (IUE), 
the United Kingdom Infrared Telescope (UKIRT), the 
James Clerk Maxwell Telescope (JCMT) in the mm wave- 
band, and the University of Michigan 26-m radio telescope 
(UMRAO) in the 4.8 - 14.5 GHz frequency range. The 
VHE 7-ray, ASCA, and EGRET observations are shown in 
Figure [HI. The VHE observations reveal the rising edge of 
a flare that developed over approximately 4 days with the 
peak flux detected on May 15 (UT) being approximately 9 
times the mean flux measured during that observing sea- 
son and approximately 1.4 times the flux of the Crab Neb- 
ula. Observations had to be halted after May 15 because 
the phase of the moon precluded further obse rvations. Ob- 
servations taken with ASCA on May 16/17 ( Takahashi et 
al. 1994 , 1996) indicated a flux approximately 20 times 
the quiescent X-ray flux of Mrk421, and the time coinci- 
dence between the two observations of unprecedented high 
states is the basis for the claim of a correlation in this cam- 
paign. Interestingly, EGRET observations taken between 
1997 May 10 and 17 did not detect the strong day-scale 
variability seen in VHE 7-rays. The average flux during 
this period was approximately twice the average flux mea- 
sured in 1994 April, so there is some evidence of a higher 
emission state, but it is not significant enough to claim 
a correlation. The observations at UV, IR, mm, and ra- 
dio wavelengths showed no evidence of variability during 
this period. Because of the offset in time of the observa- 
tions between the VHE 7-rays and the X rays, detailed 
comparisons of the variability in those wavebands are not 
possible. 

Spurred by this result, another multiwavelength cam- 
paign was organized in 1995 to better measure the mul- 
tiwavelength properties of Mrk421. This campaign re- 
vealed, for the first time, correlations between VHE 7- 



rays and X rays (Buckley et al. 1996). Observations were 
conducted between April 20 and May 5 with the Whipple 
telescope, EGRET, ASCA, the Extreme Ultraviolet Ex- 
plorer (EUVE), an optical telescope, an optical polarime- 
ter, and UMRAO. Observations with EGRET did not re- 
sult in a detection of Mrk421. The 2<r flux upper limit 
for E>100MeV is 1.2 x 10~ 7 cm -2 s _1 , somewhat be- 
low the level detected in 1994. The light curves for some 
of these observations are shown in Figure llq. The optical 
data have the contribution from the host galaxy of Mrk 421 
subtracted off. As in the 1994 multiwavelength campaign, 
Mrk 421 underwent a large amplitude flare in VHE 7-rays 
during the observations. The flare is also clearly seen in 
the ASCA and EUVE observations. There is some evi- 
dence for correlated variability in the optical flux and po- 
larization, but the statistics are not good enough to be 
confident about claiming such an association. The X rays 
and VHE 7-rays appear to vary together, limited to the 
one day resolution of the VHE observations, and the am- 
plitude of the flaring is similar, ^400% difference between 
the peak flux and that at the end of the observations. The 
EUVE and optical data (assuming it also is correlated) are 
consistent with the flare being delayed by approximately 
1 day relative to the X rays and VHE 7-rays. The am- 
plitude of the flare also decreases with decreasing energy. 
The XUV flux varies by ~200% during the observations 
and the optical flux varies by about 20%. The B-band 
percent polarization varies by nearly a factor of two in the 
observations. 

The observations of Mrk 421 in 1994 and 1995 were 
clearly undersampled, limiting the conclusions that could 
be drawn concerning correlations between wavelengths and 
emission models. Two multiwavelength campaigns orga- 
nized in 1998 attempted to improve these measurements 



18 



VHE Gamma-Ray Astronomy 



Whipple Observations 
(> 300 GeV 7-rays) 




40 



+ t + t 



(b) 



ASCA SIS0 + SIS1 

(0.7-7.5 keV X-rays) 



o 

0.4 



(c) 



/ V 






V < 



(d) 



h*' 



(e) 



*♦♦ 



EUVE DS/S 
(58-174 A XUV) 






Optical Observations 
(650 nm) 



i j + 



Optical Polarization 



T B 



49828 49830 49832 49834 49836 49838 49840 49842 49844 

Modified Julian Day 



Fig. 16. — (a) Gamma-ray, (b) X-ray, (c) extreme-UV, (d) optical, and (e) optical polarization measurements of Mrk421 taken 1995 April 
May. April 26 corresponds to MJD 49833. Figure from Buckley et al. (1996). 



through more dense observations in X rays and VHE 7- 
rays. Improvement in the VHE measurements came about 
through longer VHE exposures with individual telescopes, 
to search for hour-scale variations, and coordination of 
VHE observations between CAT, HEGRA, and Whipple. 
Thus, light curves of 12 - 16 hours could in principle be 
achieved, allowing much more detailed measurements of 
the VHE emission. These observations yielded immediate 
improvements in the measurements of flaring activity and 
correlations between VHE 7-rays and X rays. 

The first campaign, conducted in late 1998 April, 
was centered at X-ray wavelengths on observations with 
the BeppoSAX satellite and established the first hour- 
scale correlations bet ween X rays and 7-rays in a blazar 
(Maraschi et al. 1999). The lightcurve for the observations 
by BeppoSAX in three X-ray bands and Whipple above 
2 TeV is shown in Figure O. The VHE threshold here is 
higher than typical of Whipple data because observations 
were taken at a wide range of elevations, causing the en- 
ergy threshold and sensitivity to vary with time. These 
effects were corrected in the light curve through the use 
of simulations and analysis cuts to normalize the collec- 
tion area and energy threshold. Thus, the VHE lightcurve 
rate variations are intrinsic. As Figure [T| shows, a flare 
is clearly detected in X rays and TeV 7-rays on the first 
day of observations. The peaks in the lightcurves occur at 



the same time, within 1 hour, but the fall off in the X-ray 
flux is considerably slower than the TeV 7-rays. Also, the 
TeV 7-rays have a larger variability amplitude (~4-fold ra- 
tio between average and peak) than the X rays (~2-fold). 
Both the faster VHE flux decrease and the larger ampli- 
tude variability have not been seen previously in Mrk421. 
These observations provide the first clear indication that 
X rays and VHE 7-rays may not be completely correlated 
on all time-scales. 

The second multiwavelength campaign started in late 
1998 April, immediately after the observations discussed 
above, and was centered around a s even day continuous ob - 
servation of Mrk421 with ASCA ( |Takahashi et al. 1999J ). 
The light curves for the X-ray observations and VHE ob- 
servations by Whipple, CAT, and HEGRA are shown in 
Figure Il8l Had these observations been conducted as in 
1994 and 1995, with short X-ray exposures and just Whip- 
ple observations, the results would have shown nothing 
new: similar variability scale between the X rays and VHE 
7-rays and some sub-dayscale X-ray variability with am- 
plitude too low to be resolved with the Whipple observa- 
tions. Instead, the X-ray observations reveal the complete 
cycle of about 10 flares, the first time this has been done 
for Mrk421. Also, these observations seem to confirm the 
supposition of Buckley et al. (1996) that the VHE emis- 
sion from Mrk421 is primarily the result of flares, with 
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little steady emission evident. Finally, the combination of 
VHE data from these telescopes confirms the sub-day-scale 
correlations seen in the Whipple/ BeppoSAX observations. 
Detailed comparisons of the VHE 7-ray and X-ray data 
will require more sophisticated normalization of the VHE 
data (e.g., to common threshold energies) and investiga- 
tion of systematics in the measures of variability, but the 
data hold the promise of significantly advancing our study 
of VHE-cmitting 7-ray blazars. They also clearly demon- 
strate the benefits of operating multiple VHE 7-ray instal- 
lations in understanding the nature of variable sources. 

The first multi- wavelength observations of Mrk501 
which included VHE observations were conducted in 1996 
( Kataoka et al. 1999| ). The observations were conducted 
with the Whipple telescope, EGRET, ASCA, and an op- 
tical telescope. The light curve for these observations is 
shown in Figure Q9l The observations were too undersam- 
pled, or insensitive in the case of EGRET, to clearly estab- 
lish any correlations, but these observations have two very 
important results. First, follow up observations in 1996 
May established the first detection of Mrk 501 by EGRET, 
with a marginal significance of 4.0c above 100 Me V but a 
significance of 5.2cr above 500 MeV indicating a hard pho- 
ton spectrum (Kataoka et al. 1999). The claim of a 3.5a 
detection by EGRET during a small part of the multiwave- 



length campaign (the region between the dashed lines in 
Figure 19) seems speculative, given the lack of any in- 



creased activity in other wavebands during that period. 
Second, the observations established a baseline spectral en- 
ergy distribution for Mrk 501 during a relatively low emis- 
sion state which could be compared to observations of the 
high state emission in 1997 (see below for discussion of 
spectral energy distributions). 

Multiwavelength observations of Mrk 501 during its high 
emission state in 1997 revealed, for the first time, clear 
correlations between its VHE 7-ray and X-ray emission 
(Catanese et al. 1997). Observations were conducted with 



Whipple (nightly, from April 7-19), EGRET and Oriented 
Scintillation Spectrometer Experiment (OSSE) on CGRO 
(April 9 - 15), BeppoSAX (April 7, 11, 16), RXTE (twice 
nightly April 3 - 16), and the Whipple Observatory 1.2-m 
optical telescope (nightly, April 7 - 15). The optical and 
X-ray observations were serendipitously scheduled at this 
time, and the CGRO observations were a public target of 
opportunity observation initiated in response to the high 
VHE emission state. 

Figure EG shows daily flux levels for the contemporane- 
ous observations of Mrk 501. The average flux level in the 
U-band in March is also included in the figure (Fig. EGe, 
dashed line) to indicate the significant Q>10%) increase in 
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flux qctwccn March and April. An 11 day rise and fall al. 1995| ), while the Mrk501 detection had the highest 50 



in flux is evident in the VHE and X-ray wavebands, with 
peaks on April 13 and 16. The 50 - 150 keV flux detected 
by OSSE also increases between April 9 and 15, with a 
peak on April 13. The optical data may show a correlated 
rise, but the variation is small (at most 6%). Subtraction 



of the galaxy light contribution will increase the amplitude 
of this variation, but it should still remain lower than in 
X rays, given that the R-band contribut ion of the galaxy 
light is -75% (|Wurtz, Stocke k Yee 1996J) and the U-band 



contribution should be much less. EGRET observations 
indicated an excess of 1.5<r, not a significant detection. 
The ratio of the fluxes between April 13 and April 9 are 
4.2, 2.6, 2.3, and 2.1 for the VHE 7 -ray, OSSE, RXTE, 15 
- 25keV, and RXTE 2 - lOkeV emission, respectively. 

The results of this campaign show that for Mrk501, like 
Mrk 421, the VHE 7-rays and the soft X rays vary together 
and the variability in the synchrotron emission increases 
with increasing energy. However, OSS E has never detected 
Mrk 421 despite several observations ( McNaron-Brown et 



150 keV flux ever detected by OSSE from a blazar. A likely 
explanation of the OSSE detection is that the synchrotron 
emission in Mrk 501 extends to 100 keV, compared with 
the ~lkeV cutoff seen in Mrk 421. This explanatio n was 
first confirme d by the observations with BeppoSAX ( Pian 
et al. 1998 ). In addition, the day scale variations for 



Mrk 501 are larger in 7-rays than in X rays, unlike Mrk 421 . 
So, despite the similarity of Mrk 421 and Mrk 501 in some 
respects, these multi-wavelength campaigns are beginning 
to reveal differences in the two objects. 

In these short multi-wavelength campaigns, there ap- 
pears to be a correlation between the X rays and VHE 
7-rays. A natural question to ask is whether this is al- 
ways true, or only during certain situations. An attempt 
to answer this question has been made by the HEGRA 
collaboration by comparing their observations of Mrk 501 
above 500 GeV to those by the RXTE All-Sky Monitor , 
measuring 2 - 12keV photons ( Aharonian et al. 1999a ). 
A cross-correlation analysis of the daily average flux mea- 
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sured by the All-Sky Monitor with the daily average flux 
measured by HEGRA reveals a peak in the correlation 
function at At = Oil day. However, the peak in the cross- 
correlation function is only ^0.4 and the significance of the 
peak is only 2a - 3a. Whether this indicates that the X- 
ray/TeV correlation is not present is unclear because the 
ASM data have large statistical and significant systematic 
uncertainties for day-scale measurements of this relatively 
dim X-ray source (i.e., compared to the X-ray binaries the 
ASM was designed to monitor). Also, because HEGRA 
sits on the falling edge of the high energy spectrum and 
the ASM sits (for Mrk501 in 1997) on the rising edge of 
the synchrotron spectrum, it is possible that the emission 
detected by these two instruments will not be completely 
correlated, particularly for day-scale variations. Compar- 
ison of longer-term variability between the measurements 
might help resolve such issues. 

Figure pi] shows the spectral energy distributions 
(SEDs), expressed as power per logarithmic bandwidth, 
for Mrk421 and Mrk501 derived from contemporaneous 
multi-wavelength observations and an average of non- 
contemporaneous archival measurements. Both have a 
peak in the synchrotron emission at X-ray frequencies, as 
is typical of X-ray selected BL Lac objects, and a high en- 
ergy peak whose exact location is unknown but must lie in 
the 10 - 250 GeV range. Both the synchrotron and high en- 
ergy peak are similar in power output, unlike the EGRET- 
detected flat-spectrum radio quasars which can have high 
energy peaks well above the synchrotron peaks (e.g., H 



Montigny et al. 1995). Also, during flaring episodes, the 
X-ray s pectrum in both object s tends to harden signifi- 



cantly (Takahashi et al. 1996, 1999; Pian et al. 1998) 



while the VHE spectrum is not observed to change. 

The SEDs of the two sources do, however, exhibit im- 
portant differences. Most prominent among these is that 
the combination of contemporaneous RXTE and OSSE ob- 
servations of Mrk501 in 1997 clearly confirm the initial 
measurements of Pian et al. (1998) that the synchrotron 
spectrum extended well-beyond the ~lkeV typical of X- 
ray selected BL Lac objects. They also establish that the 
peak power output of the synchrotron emission occurs at 
~100keV. This is in contrast to the 1996 observations of 
Mrk501 reported by Kataoka et al. (1999), where the 
synchrotron power peak is at ~2keV. For Mrk421, the 
X-ray spectral peak does shift to higher energies during 
flaring activity, but the changes are much smaller than 
in Mrk501, and the peak was never observed to extend 
beyond ~lkeV. This peak is followed by a sharp cutoff 
which produces a deficit in the OSSE range, preventing 
the detection of Mrk421 by this instrument. 

Whether the shift in the location of the synchrotron 
peak for Mrk501 is also accompanied by a shift in the on- 
set of the 7-ray emission to higher energies is not clear. 
Any increase in the MeV-GeV flux in 1997 was not as 
great as at TeV energies, or EGRET would have easily 
detected Mrk501. But, because the sensitivity of EGRET 
in 1997 was substantially poorer than in 1996, a shift in 
the onset of the spectrum or a flux variation that increases 
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Catanese (1999). 



with increasing energy could explain the non-detection by 
EGRET. The fact that the VHE spectrum of Mrk 501 is 
harder than that of Mrk 421 below ~1 TeV (see Figure ||) 
may also indicate that the high energy peak of Mrk 501 is 
shifted to slightly higher energies. However, it could also 
simply indicate a slower fall off in the progenitor particle 
spectrum above the peak power output. 

A second difference in the spectral energy distributions 
for these two objects is that the power output for Mrk 501 
in the VHE range can be considerably less than in X rays 
when it is in a low emission state. In contrast, Mrk 421 
seems to maintain a similar output at X-ray and 7-ray en- 
ergies. These differences are illustrated in Table a which 
gives the ratio of contemporaneously measured fluxes for 
X rays and 7-rays for these two objects. Whether the dif- 
ference in power output for Mrk 501 reflects a change of 
the energy at which the peak in the high energy spectrum 
occurs or something related to the flaring process is not 
clear, due to the poor spectral measurements in the low 
emission states and the lack of coverage of the peak region 
of the spectrum. However, the lack of spectral variability 



in 7-rays argues against a significant short-term shift of 
the 7-ray spectral peak. 

4.1.3. Implications of the Very High Energy Observations 

The general properties of the detected extragalactic 
sources of VHE 7-rays are listed in Table g. The three 
objects detected by the Whipple Collaboration exhibit 
some interesting commonalities. They are the three clos- 
est known BL Lac objects with declination > 0° so their 
7-ray fluxes are the least attenuated from interaction with 
background IR radiation. The >100MeV fluxes are near 
(Mrk 421, Mrk 501) or below (1ES 2344+514) the EGRET 
sensitivity limit, meaning that the 7-ray power output 
does not peak in that energy range as it does for many o f 
the EGRET-detected AGNs Qvon Montigny et al. 1995| ). 
Thus, VHE observations already augment the catalog of 
7-ray sources compiled by space-borne telescopes. Fi- 
nally, all three of the Whipple-detected BL Lac objects 
are X-ray selected BL Lac objects (XBLs). The exten- 
sion of the synchrotron spectra to X-ray energies in XBLs 
implies that they produce high energy electrons, making 
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Fig. 21. — Left: The spectral energy distribution of Mrk421 from contemporaneous and archival observations. Dates of the observations 
are indicated in the figure. Figure from Buckley et al. (1997). Right: The spectral energy distribution of Mrk501 from contemporaneous and 
archival observations. Dates of the observations are indicated in the figure. Data in the figure come from Kataoka et al. (1999), Catanese 
et al. (1997), and Catanese (1999) and references therein. The curves in the figure are meant to guide the eye to the contemporaneously 
measured points, and do not indicate model fits to the data nor are they an attempt to elucidate the spectral energy distribution of Mrk 501 
during these observations. In both figures, the archival measurements are approximate averages of the data in the literature. 

Table 5 
Spectral energy ratios for VHE sources 



Source 


Date 


pa 
-"^kcV 


-"rookcv 


-*noo McV 


Rcf. 


Mrk 421 


1995 
1996 
1996 
1997 


2.2±0.5 

0.81±0.06 

4.4±2.1 

1.2±0.3 


1.9±0.5 


<0.38 

3.1±3.4 
<0.36 




Buckley et al. 1996 






Buckley et al. 1997 




Mrk 501 


Kataoka et al. 1999 






Catanese et al. 1997; Catanese 199E 



i/F I/ (2kcV)/z/F„(350GeV) 
b i/F iy (100keV)/^F iy (350GeV) 
c z/F^(100MeV)/i/F I/ (350GeV) 



them good candidates for VHE emission if the VHE 7- 
rays are produced via inverse Compton (IC) scattering of 
these same electrons. EGRET's tendency to detec t more 
radio- sele cted BL Lac objects (RBLs) than XBLs ( Lin et 
al. l j)97] ) also supports this tenet because RBLs would 
be exp ected to have spect r a which peak in the MeV -GeV 
range flSikora et al. 1994J; JMarschcr fc Travis 19961). BL 



Lac objects in general have been suggested as better can- 
didates for VHE emission than other blazars because the 
absence of optical emission lines in BL Lac objects may 
indicat e less VHE-absorbing radiat ion near the emission 
region ( Dcrmcr fc Schlickciser 1994 ). 

The other two objects detected at VHE energies, 
PKS 2155-304 and 3C66A, are similar in some respects 
to the Whipple sources. PKS 2155-304 is an XBL, so it 
fits the IC paradigm for VHE sources. However, 3C 66A is 
classified as an RBL, suggesting that protons produce the 
7-rays because in IC models, RBLs would not have high 
enough energy electrons to produce TeV emission. In ad- 
dition, both PKS 2155-304 and 3C 66A are at much higher 
redshifts than the Whipple sources, implying quite low IR 
backgrounds. Thus, confirmation of these detections, just 



as for 1ES 2344+514, is essential. 

VHE observations have already significantly affected our 
understanding of BL Lac objects. For example, the rapid 
variability indicates either very low accretion rates and 
photon densities near the nucleus (Celotti, Fabian & Rees 
1998) or, conversely, requires the 7-ray emission region to 
be located rel atively far from the nucleus to escape the 
photon fields ( Prothcroc fc Bicrmann 1997 ). Also, the ob- 



servations have helped resolve the nature of the differences 
between RBLs and XBLs. Based on their smaller numbers 
and higher luminosities, Maraschi et al. (1986) proposed 
that RBLs were the same as XBLs but with jets aligned 
more closely with our line of sight. However, the rapid 
variability and TeV extent of the XBL emission point to 
the differences between the two sub-classes being more fun- 
damental, as originally proposed by Padovani & Giommi 
(1994): the XBLs have higher maximum electron energies 
and lower intrinsic luminosities. 

Simultaneous measurements of the synchrotron and 
VHE 7-ray spectra also constrain the magnetic field 
strength (B) and Doppler factor (S) of the jet. If the 
correlation between the VHE 7-rays and optical/UV pho- 
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Table 6 
Properties of the VHE BL Lac objects 







EGRET fmx a 


Average flux 




^x a 


^R a 


Object 


z 


(EMOOMeV) 


(E>300GeV) 




(2keV) 


(5 GHz) 






(10- 7 cm- 2 s- f ) 


(lO-^cm^s- 1 ) 


M v a 


(MJy) 


(mJy) 


Mrk421 


0.031 


1.4±0.2 


40 


14.4 


3.9 


720 


Mrk 501 


0.034 


3.2±1.3 


>8.1 


14.4 


3.7 


1370 


1ES2344+514 C 


0.044 


<0.7 


^8.2 


15.5 


1.1 


220 


PKS 2155-304° 


0.116 


3.2±0.8 


42 


13.5 


5.7 


310 


3C66A C 


0.444 


2.0±0.3 


30 6 


15.5 


0.6 


806 



a Radio, optical, and X-ray data from Pcrlman et al. (1996). EGRET data from D.J. Thomp- 

son (priv. comm.), Mukherjee et al. (1997), and Kataoka et al. (1999). 

b >lTeV flux value. 

c Unconfirmed as a VHE source. 



tons observed in 1995 from Mrk 421 indicates both sets 
of photons are produced in the same region of the jet, 
S ;> 5 is required for the VHE photons to escape signifi- 
cant pair- production losses (Buckley et al. 1996). If the 
SSC mechanism produces the VHE 7-rays, 6 — 15 — 4 
and B = 0.03 - 0.9G for Mrk 421 flBuckley et al. 199 
Tavecchio, Maraschi & Ghisellini 1998; |Catancsc 1999 1 
and(5 « 1.5-20 land ff = 0.08-0. 2G for Mrk 501 ([jamuc , 
son el j al. 1998; ; |Tavecchio et al. 1998J ; Hillas 1999| )~~To 



match the variability time-scales of the correlated emis 
sion, proton models which utilize synchrotron cooling as 
the primary means for proton energy lo sses require mag- 
netic fields of B = 30 - 90G for 5 « 10 ( |Mannhcim 1993 ; 
Mannheim 1998J ; |Bucklcy 1998J ). The Mrk 421 values of 
6 and B are extreme for blazars, but they are still within 
allowable ranges and are consistent with the extreme vari- 
ability of Mrk 421. 

In addition, the VHE observations have constrained the 
types of models that are likely to produce the 7-ray emis- 
sion. For instance, the correlation of the X-ray and the 
VHE flares is consistent with IC models where the same 
population of electrons radiate the X rays and 7-rays. The 
absence of fla ring at EGRET ener gies may also follow in 
this context ( Macomb et al. 1995] ) because the lower en- 
ergy electrons which produce the 7-rays in the EGRET 
range radiate away their energy more slowly than the 
higher energy electrons which produce the VHE emission. 
The MeV-GeV emission could then be the superposition 
of many flare events and would therefore show little or no 
short-term variation. 

In the mechanism of Sikora et al. (1994), which produces 
7-rays through the Comptonization of external photons, 
the external photons must have energies <0.1eV (in the 
IR band) to avoid significant attenuation of the VHE 7- 
rays by pair production. Sikora et al. (1994) point out that 
there is little direct observational evidence of such an IR 
component in BL Lac objects, but the existence of such a 
field has been pre dicted as a product of accretion in AGNs 
( |Rees et al. 1982; ). 

Models which produce the 7-ray emission from proton 
progenitors through e + e~ cascades originating close to the 
base of the AGN jet have great difficulty explaining the 
TeV emission observed in Mrk 421 because the high den- 



sities of unbeamed photons near the nucleus, such as the 
accretion disk or the broad line region, required to initi- 
ate the cascades cause high pair opacities to TeV 7-rays 
(Coppi, Kartje & Konigl 1993). Such models also predict 
that the radius at which the optical depth for 7-7 pair 
production drops below unity increases with increasing 7- 
ray energy (Blandford & Levinson 1995) and therefore the 
VHE 7-rays should vary either later or more slowly than 
the MeV-GeV 7-rays. This is in contradiction to the ob- 



servations of^lrjc421Jnj50th 1994 ( [Macomb et al. 1995[ ) 
and 1995 (|Bucklcy et al. 1996J ). 



4.2. Extragalactic Background Light 

In traversing intergalactic distances, 7-rays may be ab- 
sorbed by photon-photon pair production (7 + 7 — > e + + 
e~) on background photon fields if the center of mass en- 
ergy of the photon-ph oton system exceeds tw ice the rest 
energy of the electron (Gould & Schreder 1967). The cross- 
section for this process peaks when 



E 7 e(l 



cos 



9) ~ 2(m e c 2 ) 2 = 0.52(MeV) 2 



(1) 



where E 7 is the energy of the 7-ray, e is the energy of the 
low energy photon, 9 is the collision angle between the 
two photons, m e is the mass of the electron, and c is the 
speed of light in vacuum. Thus, for photons of energy near 
1 TeV, head-on collisions with photons of ^0.5 eV have the 
highest cross-section, though a broad range of optical-to- 
IR wavelengths can be important absorbers because the 
cross-section for pair production is rather broad in energy 
and spectral features in the extragalactic background den- 
sity can make certain wavebands more important than the 
cross-section alone would indicate. 

The presence of extragalactic background light (EBL) 
limits the distance to which VHE 7-ray telescopes can de- 
tect sources. This has been put forth as an explanation 
of the lack o f detection of many of the EGRET- detected 
AGNs (e.g., Stecker, de Jager fc Salamon 1992), as dis- 



cussed above. The difficulty in understanding the effect of 
the EBL on the opacity of the universe to VHE 7-rays is 
that not much is known about the spectrum of the EBL 
at present, nor how it developed over time. Star forma- 
tion i s expected to be a major contributor t o the EBL 
(e.g., Madau et al. 1996: Primack et al. 199E), with star 
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formation contributing mainly at short wavelengths (1 - 
15 /im) and dust absorption and re-emission contributing 
at longer wavelengths (15 - 50 /mi) . So, measurements 
of the EBL spectrum can serve as important trac ers of 
the h istory of the formation of stars and galaxies ( Dwek 



et al. 1998| ). Other, more exotic processes, such as pre 



galactic star formation and some dark matter candidates, 
might also contribute distin ctive features to the EBL (e.g., 
Bond, Carr fc Hogan 1986 , 1991). Thus, measurements of 
the EBL have the potential to provide a wealth of infor- 
mation about several important topics in astrophysics. 

Experiments that attempt to measure the EBL by di- 
rectly detecting optical-IR photons, such as the Diffuse 
Infrared Background Experiment (DIRBE) on the Cosmic 
Background Explorer ( COBE) , are plagued by foreground 
sources of IR radiation. Emitted and scattered light from 
interplanetary dust, emission from unresolved stellar com- 
ponents in the Galaxy, and dust emission from the inter- 
stellar medium are all significantly more intense than the 
EBL and must be carefully modelled and subtracted to de- 
rive estimates of the EBL. Currently, EBL detections ar e 
available only at 140 /jm and 24 /im (Hauser et al. 1998 ). 
Tentative detec tions at 3.5 ^m (Dwek fc Arendt 1998 ) and 
400 - 1000 fjja (Puget et al. 1996) have also been reported. 

Because VHE 7-rays are attenuated most by optical-IR 
photons, measurements of the spectra of AGNs provide 
an indirect means of investigating the EBL that is not af- 
fected by local sources of IR radiation (|Gould fc Schreder 



1967f| |Stecker, de Jager fc Salamon 1992J ). The signs of 



EBL absorption can be cu toffs, but also simple alterations 
of the spectral index (e.g., Stecker 1999), depending on the 
spectral shape of the EBL and the distance to the source. 
Like direct measurements of the EBL, this technique has 
difficulties to overcome. For instance, it requires some 
knowledge of or assumptions about the intrinsic spectrum 
and flux normalization of the AGNs or the EBL. Also, the 
AGNs themselves produce dense radiation fields which can 
absorb VHE 7-rays at the source and thereby mimic the 
effects of the intergalactic EBL attenuation. 

Despite these difficulties, the accurate measurement of 
VHE spectra with no obvious spectral cut-offs from just 
the two confirmed VHE-emitting AGNs, Mrk421 and 
Mrk501 (see § [l.l.l| ), has permitted stringent limits to 
be set on the density of the EBL over a wide range of 
wavelengths. These limits have been derived from two 
approaches: (1) assuming a limit to the hardness of the 
intrinsic spectrum of the AGNs and deriving limi ts which 
assume v ery little about the EBL spect rum (e.g., Billcr et 



al. 1£ 9S ; Stanev fc Franccschini 1998|) and (2) assuming 



some shape for the EBL spectrum, based on theoretical 
or phenomenological modelling of the EBL, and adjust- 
ing the normalization of the EBL density to match the 
mea sured VHE spectra (e.g., ie Jager, Stecker fc Salamor] 



1994 | [Btanev fc Franceschini 1998 ). The latter can be more 
stringent, but are necessarily more model-dependent. The 
limits from these indirect methods and from the direct 
measurements of EBL photons are summarized in Fig- 
ure E2. At some wavelengths, the TeV limits represent 
a 50-fold improvement over the limits from DIRBE. These 
limits are currently well above the pre dicted density for the 
EBL from normal galaxy formation (Madau et al. 1996 



EBL (e.g., pillcr et al. 1998| ). They also show that EBL 
attenuation alone cannot explain the lack of detection of 
EGRET sources with nearby redshifts at VHE energies, as 
the optical depth for pair-production does not reach 1 for 
the stringent limits of Biller et al. (1998) until beyond a 
redshift of z = 0.1 (see Figure £3|). With the detection or 
more AGNs, particularly at higher redshift, and improve- 
ments in our understanding of the emission and absorption 
processes in AGNs, VHE measurements have the poten- 
tial to set very restrictive limits on the EBL density, and 
perhaps eventually detect it. 

4.3. Gamma-Ray Bursts 

Although the 7-ray burst phenomenon is usually asso- 
ciated with energies of 100 keV to 1 MeV (hard X rays 
to low energy 7-rays) results from EGRET show that 
there is a component at high energies and thus the phe- 
nomenon has the potential to be observed in the TeV 
range. The power spectrum certainly peaks in the lower 
energy ranges but the observations at high energies really 
provide the strongest constraints on the emission models 
and may ultimately expose the underlying emission mech- 
anism. The detection of a single photon of energy 18 GeV 
from GRB 970217, 1.5 hours after the onset of the burst 



( Hurley et al. 1994 ) , has opened th e possibility of delayed 
emission of GcV -TeV photons (e.g., Totani 1998 ; Bottchci 
fc Dermer 1998| ) . Althoug h ACITs have, to date, onl y pre- 
sented upper limits (e.g., |Connaughton et al. 1997 ), it is 
possible that in the near future wide field air shower de- 
tectors like MILAGRO or the Tibet Array might detect a 
prompt VHE emission component and rapid slew ACITs 
might do the same with greater sensitivity for any delayed 



Primack et al. 1999| ), but they have provided constraints 
on a variety of more exotic mechanisms for sources of the 



emission. 

5. THE FUTURE OF VHE ASTRONOMY 

It is clear that to fully exploit the potential of ground- 
based 7-ray astronomy the detection techniques must be 
improved. This will happen by extending the energy cover- 
age of the technique (with good energy resolution) and by 
increasing its flux sensitivity (improved angular resolution 
and increased background rejection). Ideally one would 
like to do both but in practice there must be trade-offs. 
Reduced energy threshold can be achieved by the use of 
larger, but cruder, mirrors and this approach is currently 
being expl oited using existing arrays of sola r hcliostats: 
STACEE dChantcll et al 199 j), CELEST E dQuebert et 
al. 1995| ) and Solar-2 ( Turner et al. 1999|). A German- 
Spanish project (MAGIC) ( [Barrio et al. 1998| ) to build 
a 17-m aperture telescope has also been approved. These 
projects will achieve thresholds as low as 20-30 GeV where 
they will effectively fill the current gap in the 7-ray spec- 
trum from 20 to 200 GeV. Ultimately, this g ap will be cov- 
ered by GLAST ( |Gchrcls fc Michclson 1999J ). Extension to 
higher energies (>10 TeV) can be achieved by atmospheric 
Cherenkov telescopes working at large zenith angles and 
by particle arrays on very high mountains. An interesting 
telescope that has just come on line and will complement 
these techniques is the MILAGRO water Cherenkov de- 
tector in New Mexico which will operate 24 hours a day 
with a large field of view an d will have good se nsitivity to 
7-ray bursts and transients ( [Binnis et al. 1995 ). 

VERITAS, with seven 10-m telescopes arranged in a 
hexagonal pattern with 80 m spacing (Figure 04) , will aim 
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for the middle ground between those techniques listed 
above, with its primary objective being high sensitivity 
observations in the 100 GeV to lOTeV range ( Weekes et 
al. 19w9f). It will be located in southern Arizona and will 



be a l pgical progression from the Whipple telescope. It is 
hoped to begin construction in 1999 and to complete the 
array by 2004. 

The German- French-Italian experiment HESS, initially 
four and e ventually perha ps sixteen 12m class telescopes in 
Namibia ( Hofmann 1997 ), and the Japanese Super CAN- 
GAROO array, with four 10-m telescopes in Australia, (T. 
Kifune, private communication) will have similar objec- 
tives. In each case, the arrays will exploit the high sensitiv- 



ity of \CITs and the high selectivity of the array approach. 
The projected sensitivities of MAGIC, HESS, SuperCAN- 
GAROO and VERITAS are somewhat similar and we re- 
fer to them collectively as Next Generation Gamma-Ray 
Telescopes (NGGRTs). The relative flux sensitivities for 
existing and planned 7-ray telescopes as a function of en- 
ergy are shown in Figure 25 , where the sensitivities of the 
wide field detectors are for one year and the atmospheric 
Cherenkov telescopes are for 50 hours. In all cases, a 5 a 
point source detection is required. 

It is apparent from this figure that, on the low energy 



side (<lTeV), the NGGRTs will complement the GLAST 
mission and will overlap with the solar arrays. At the 
highest energies to which they are sensitive, NGGRTs will 
overlap with the Tibet Air Shower Array (Amenomori et 
al. 1997). They will cover the same energy range as MILA- 



GRO but with greater flux sensitivity. The wide field cov- 
erage of MILAGRO will permit the detection of transient 
sources which, once detected, can be studied in more de- 
tail by the northern NGGRTs. These same telescopes will 
complement the coverage of neu trino sources to be discov- 
ered by AMANDA/ICE CUBE ( |Halzcn 199&J ) at the South 
Pole. Finally, if the sources of ultra-high e nergy cosmic 



rays are loc alized to a fe w degrees by H iRes (Abu-Zayyad 



et al. 1997) and Auger (Boratav 1997), the NGGRTs will 



be powerful instruments for their further localization and 
identification. 

The recent successes in VHE 7-ray astronomy en- 
sure that in the inevitable interval between the death of 
EGRET and the launch of the next generation 7-ray space 
telescope, there will be ongoing activity in GeV-TeV 7-ray 
astronomy. Observations by GLAST and the NGGRTs 
in this energy region will make important contributions 
to our understanding of AGNs, supernova remnants, and 
pulsar and 7-ray burst studies. Although the number of 



Catanese & Weekes 



27 



10° 



m 
a 



10' 



10' 




10' 



10" 1 10° 

Energy (TeV) 



10 1 



Fig. 23. — Lower limit of the redshift at which photons of the energies shown face an optical depth of 1 due to pair-production with EBa.L 
Limits are derived from the upper limits on the density of the EBL of Biller et al. (1998). Figure from Biller et al. (1998). 




Fig. 24. — The proposed arrangement of telescopes in VERITAS. Figure from Weekes et al. (1999). 



TeV sources detected so far is small, the new and varied 
phenomena observed indicate that VHE 7-ray astronomy 
is not merely an extension of McV-GcV 7-ray astronomy 



but is a discipline in its own right. With the advent of 
new telescopes, the catalog of VHE 7-ray sources will dra- 
matically expand with detailed time histories and accurate 
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Fig. 25. 
CFLESTE/STACFF flrjucbcrt et al. 1995| ; 

( [Thompson et al. 1993[ ), and MLLAURU ( fcinms et al. 19951 ).' The sensitivity ot 

e.g., hybrid PMTs, not assumed in the other experiments. EGRET, GLAST and MILAGRO are wide field instruments and therefore ideally 
suited for all sky surveys. The turn-up in the VERITAS sensitivity at higher energies is primarily caused by the requirement that the signal 
contain at least 10 photons. 
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